lag <- function(x, period, Nlag, N){
       if (Nlag >= 0){
              tmp <- rep(NA, N)
	      tmp[period] <- x
	      x <- tmp       
	      N <- length(x)	   
              out <- c(rep(NA, Nlag), x[1:(N-Nlag)])
       } else {
              tmp <- rep(NA, N)
	      tmp[period] <- x
	      x <- tmp       
	      N <- length(x)
	      Nlag <- abs(Nlag)	   
              out <- c(x[(1 + Nlag):N], rep(NA, Nlag))
       }       
       return(out)
}
